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Abstract 

We test an isospectral potential from harmonic oscillator 
simulating H-bond interation in DNA macromolecules. 


1 Introduction 

In the context of Supersymmetric Quantum Mechanics (SQM), several new potentials have 
been generated from shape invariant potentials[l] . In particular, with reference to the 
harmonic oscillator[2-4]. However, to date, none of these new potentials have been applied 
to real physical systems. 

Peyrard and Bishop[5] have introduced a theoretical treatment of DNA. In this model, 
starting from the classical Hamiltonian and using the transfer integral operator, the partition 
functions follow from the eigenvalues of a Schrodinger-type equation. 

We intend, in this communication, to demonstrate the feasability of using a new potential 

i 

generated from the harmonic oscillator in the Schrodinger-type equation to describe H- 
bonds in DNA. 

2 SQM 

In SQM the nilpotent operators Q and Q + satisfy the algebra[6] 

{Q.Q + } = H SS ;Q 2 = Q +2 = 0 (1) 
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which can be realized in a usual matrix form. In this case, the supersymmetric Halmitonian 
H ss consists of two partners H + and H. which can be factorized in H + = a + a‘ and H. = a' a + 
(where a + 5^^^J^^Jbosonic operators), and they have the same spectra except for the 
ground state, where only H+ has an eigenvalue equal to zero E+ = 0. The eigenfunctions of 
H + and H. are related to one another by: 4*+ = a + 4*. and 4*. = a" 4 / + , where the bosonic 
operators are defined by: 

a ± = {+j- + W(x)}. (2) 

dx 

Redefining the operators a 1 , it is possible to obtain a family of Hamiltonians as shown 
by Sukumar[7] , or the isoespectral potentials[2-4] . In this second case defining new 
operators: 


A ± =T~ + F(x) (3) 

dx 

and imposing H. = A" A + yields a general form to get F(x). The new Hamiltonian is written 
as 


H + = A“A + - 2-^— F(x) (4) 

dx 

and the eigenfunctions of the new Hamiltonian are related to the original Hamiltonian ones 
. The missing ground state is obtained from A' 4\ 0 = 0. 

As an example, consider the original potential to be the harmonic oscillador, V(x)=x . It 
then follows[ 1-4] : 

( 5 ) 

( 6 ) 

3 DNA Model 


~ d 2 7 d exp(-x 2 ) 

H. = - — j+ x 2 - 1-2 
+ dx 2 


T + J*exp(-z 2 )dz 


7 fx exp(-z )dz 

%,0 = exp{-x 2 + I £ ) . 

["+ J j exp(-y 2 )dy 
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Bishop & Peyrard's DNA model includes two degrees of freedom (u n ,v n ) corresponding to 
displacements of the bases from their equilibrium positions along the direction of the H- 
bonds connecting them. A harmonic coupling due to the stacking interactions between 
neighboring bases is also assumed with the same coupling constant (k) for the two strands. 
Each base has the mass m. 

The Hamiltonian for this model is 


h = + |(x„ - x n _,) 2 + ^ + |(y«i - y n -i ) 2 + v (y n )} » ( ? ) 

where x n = (u n + v n )Ui , y n = (u n - v n )Ui , p n = nu„, q n = m y n and V(y n ) is the potential 
for the hydrogen bonds. Using the transfer integral operators to solve configuration integral 
energy the classical partition function can be obtained. 

The temperature of denaturation of DNA can be monitored through the mean stretching 
<y n >of the H-bond. Considering the limit of large N, only the ground state will be 
important, and 

< y >=< d> 0 (y)|y|<D 0 (y) >= jo 0 2 (y)ydy , (8) 

where <J> 0 (y) is the ground state eigenfunction of a Schrodinger-type equation: 


{- 


1 d 2 
2p 2 k dy 2 


+ V(y)}O i (y) = {s i + 


1 2n 

2pV W5,) ’ 


( 9 ) 


It is important to note the mean value < y > depends on the form of the eigenfunction 

d> 0 (y)- 


4 H-bond Potential 

The Morse potential is usually used to simulate the H-bonds in DNA. However, we suggest 
here that the potential generated from the harmonic oscillator using the superalgebra can 
also be used. This potential has the form: 
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which can be realized in a usual matrix form. In this case, the supersymmetric Halmitonian 
H ss consists of two partners H + and H. which can be factorized in H+ = a + a’ and H. = a' a + 
(where a + and a' are bosonic operators), and they have the same spectra except for the 
ground state, where only H + has an eigenvalue equal to zero E+ = 0. The eigenfunctions of 
H + and H. are related to one another by: V P+ = a + T, and T_ = a‘ 4*+ , where the bosonic 
operators are defined by: 

a± = { +£ + W( x)}. (2) 

Redefining the operators a* , it is possible to obtain a family of Hamiltonians as shown 
by Sukumar[7] , or the isoespectral potentials[2-4] . In this second case defining new 
operators: 

A±=T £ + F(X) (3) 

and imposing H. = A' A + yields a general form to get F(x). The new Hamiltonian is written 
as 


H + =A-A + -2^F(x) (4) 

and the eigenfunctions of the new Hamiltonian are related to the original Hamiltonian ones 
. The missing ground state is obtained from A" = 0. 

i , 

As an example, consider the original potential to be the harmonic oscillador, V(x)=x . It 
then follows [1-4] : 


fL = — — y+x 2 -l-2-^ eXP( ~ X) 

dx 


‘ +_ ' "dx 


r+ 


_ s 

.*exp(-z 2 )dz 


( 5 ) 


%,-exp' ■ >. 

* r+ J ) exp(-y 2 )dy 


( 6 ) 
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